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Mid-Frequency  Sonar  Interactions  with  Beaked  Whales: 

Final  Report 
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Phone:  (508)  289-2445  FAX:  (508)  457-2 1 94  E-mail:  k footer  w hoi.edu 

Award  Number:  N000 1407 10992 


EXECUTIVE  SUMMARY 

Naval  exercises  involving  use  of  military  sonars  operating  in  the  mid-frequency  (MF)  band,  1-10 
kHz,  have  been  implicated  in  the  strandings  of  marine  mammals,  especially  beaked  whales. 
Understanding  the  cause  or  causes  of  the  greater  sensitivity  of  beaked  whales  to  MF  sonar  has 
stimulated  much  speculation  and  thought.  The  suspicion  that  there  is  a  deep  connection  to 
anatomy  instigated  the  project  summarized  here,  namely  “Mid-frequency  sonar  interactions  with 
beaked  whales.”  The  ambition  of  this  was  to  develop  a  tool  that  would  enable  researchers  to 
predict  sonar-induced  acoustic  fields  inside  beaked  whales  as  well  as  inside  other  marine 
mammals.  The  tool,  called  the  Virtual  Beaked  Whale,  was  to  be  an  interactive  online  modeling 
and  visualization  system.  Formidable  resources  were  contributed  to  development  of  this  system 
both  by  ONR  and  by  the  Woods  Hole  Oceanographic  Institution,  but  it  was  not  completed. 
However,  the  scope  and  components  of  such  a  system  were  outlined  in  detail,  and  some  of  these 
components  were  also  developed.  This  included  specification  of  a  finite-element  method  (FEM) 
to  solve  the  full-diffraction  wave  equation  for  acoustic  interactions  with  a  complicated, 
heterogeneous  body;  modeling  the  three-dimensional  anatomy  of  exemplary  marine  mammal 
specimens;  and  assigning  physical  properties  to  the  diverse,  constituent  tissues.  Beaked  whales 
are  large:  a  mature  Cuvier’s  beaked  whale  ( Ziphius  cavirostris)  can  be  8.5  m  in  length.  This 
presents  two  difficulties.  Firstly,  detailed  anatomical  data  cannot  be  obtained  on  a  whole 
specimen  by  computed  tomography  (CT)  scanning.  Thus,  a  synthetic  model  of  a  beaked  whale 
was  developed  based  on  CT  scans  of  (i)  the  head  of  a  freshly  dead  Cuvier’s  beaked  whale,  and 
(ii)  the  body  of  a  harbour  porpoise  ( Phocoena  phocoena).  Secondly,  given  MF  wav  elengths  of 
order  0.15-1 .5  m,  the  total  number  of  equations  to  be  solved  is  of  order  several  tens  of  millions  at 
the  upper  end  of  this  band.  This  is  a  daunting  number.  Consequently,  particular  attention  was 
paid  to  validation  of  the  FEM  model.  For  empirical  validation,  an  experiment  was  performed  at 
the  Naval  Surface  Warfare  Center,  Cardcrock  Division,  in  which  an  instrumented  common 
dolphin  ( Delphimis  delphis)  cadaver  w'as  ensonified  and  internal  fields  measured  in  a  number  of 
tissues.  In  the  course  of  this  work  the  discovery  was  made  that  tourmaline  sensors,  which  are 
ordinarily  used  to  measure  acoustic  pressures  at  explosive  shock  levels,  can  also  be  used  to 
measure  weak,  sub-shock  pressures.  For  theoretical  validation,  provision  was  made  for 
benchmarking  through  comparison  with  analytic  solutions  for  internal  acoustic  and  particle 
displacement  fields  due  to  ensonification  of  lossy  fluid  and  solid  elastic  spheres  by  plane 
harmonic  waves.  Ultimately,  the  validated  model  was  to  serve  as  the  heart  of  the  Virtual  Beaked 
Whale.  Preliminary  work  was  done  on  developing  the  user  interface. 

The  final  project  report  mostly  follows  the  format  of  annual  NOPP  and  ONR  project  reports. 
After  a  formal  introduction,  the  long-term  goals,  objectives,  and  approach  are  presented 
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succinctly.  The  completed  work  is  then  summarized  by  task  and  year.  The  principal 
accomplishments  of  the  project  are  listed  in  the  conclusions  section  together  with  comments  on 
potential  future  work  based  on  present  findings,  ultimately  to  understand  physical  effects  of 
sonar  interactions  with  beaked  whales. 

INTRODUCTION 

Mid-frequency  (MF)  military  sonars,  with  operating  frequencies  in  the  1-10  kHz  band,  have 
been  implicated  in  a  number  of  strandings  of  beaked  w  hales.  These  strandings  have  included 
documented  events  in  Greece  in  1996  (Frantzis  1998,  D'Spain  et  al.  2006),  in  the  Bahamas  and 
Madeira  in  2000  (NOAA  2001,  Cox  et  al.  2006),  and  in  the  Canary  Islands  in  2002  (Proceed. 

ECS  2004).  Cuvier’s  beaked  whale  (Ziphins  cavirostris)  w;as  especially  affected.  Some  later 
mass  stranding  events  were  also  coincident  w  ith  the  use  of  MF  sonar,  as  in  the  Haro  Strait  near 
Washington  State  in  2003  (Norman  et  al.  2004,  National  Marine  Fisheries  Service  2005)  and  off 
the  coasts  of  Hawaii  in  2004  (Southall  et  al.  2006)  and  North  Carolina  in  2005  (Flohn  et  al. 

2006).  These  later  events  involved  other  species,  but  without  the  kinds  of  injuries  reported  in 
beaked  whales  (Freitas  2004,  Ketten  2005,  Fernandez  et  al.  2005).  This  historv  indicates  a 
particular  sensitivity  of  beaked  whales  to  MF  sonar.  It  has  motivated  a  number  of  studies, 
including  the  present  one. 

Notable  hy  potheses  that  have  been  examined  include  those  of  acoustic  resonance  in  the  lungs 
and  formation  of  gas  bubbles  in  potentially  super-saturated  deep-diving  mammals.  The  first  of 
the  hypotheses  was  eventually  discounted  (NOAA  2002). 

Research  since  the  mid-1990s  has  focused  on  sound  source  characteristics  and  propagation  paths, 
location  and  abundance  of  marine  mammal  populations,  criteria  and  thresholds  for  physical  and 
behavioral  effects  (National  Research  Council  2000,  2003,  2005), and  morphology-based 
modeling  (Aroyan  2001,  Krysl  et  al.  2006,  Cranford  et  al.  2008a).  The  models  are  realizations  of 
the  finite-element  method  (Burnett  1987). 

It  was  the  ambition  of  the  project  summarized  here  to  develop  a  tool  that  w  ould  enable  the 
community  to  investigate  the  physical  effects  of  sonar-induced  sound  on  beaked  whales  and 
other  marine  mammals.  This  tool  would  be  an  interactive  online  modeling  and  visualization 
system.  It  would  enable  researchers  to  predict  sonar-induced  acoustic  fields  inside  beaked 
whales  as  well  as  inside  other  marine  mammals. 

In  the  event,  development  of  the  tool,  called  the  Virtual  Beaked  Whale,  w  as  not  completed. 
However,  several  other  goals  were  achieved,  including  synthesis  of  a  morphometric  model  of  a 
beaked  whale.  This  and  work  on  other  tasks,  as  well  as  achievements,  are  summarized  below. 

The  basic  form  of  this  report  follows  that  of  annual  project  reports.  Long-term  goals  and 
scientific  and  technological  objectives  are  stated  succinctly,  and  the  intended  approach  is 
summarized.  Completed  work  is  then  described  for  each  task  by  year,  essentially  quoting  from 
the  previous  annual  reports  but  updating  the  information  too.  Main  deliverables  of  the  project 
are  listed  in  the  conclusions,  which  also  point  to  potential  future  work  based  on  the  present 
findings. 


LONG-TERM  GOALS 

The  top-level  goal  of  this  project  was  to  build  an  interactive  online  modeling  and  visualization 
system,  called  the  Virtual  Beaked  Whale,  to  enable  users  to  predict  mid-frequcncy  sonar-induced 
acoustic  fields  inside  beaked  whales  and  other  marine  mammals.  Another  high-level  goal  was  to 
acquire  new  high-resolution  morphometric  and  physical-property  data  on  beaked  whales  for  use 
in  the  model.  It  was  hoped  that  the  availability  of  such  a  system  together  with  high-quality  data 
would  give  researchers  insight  into  the  nature  of  sonar  interactions  with  beaked  whales, 
ultimately  to  introduce  objectivity  into  a  public  discussion  that  has  been  hampered  by  lack  of  a 
scientific  approach.  It  was  hoped  further  that  the  tool  would  prove  useful  in  evaluating  alternate 
sonar  transmit  signals  that  retain  the  required  information  content  but  with  substantially  reduced 
physical  effects  in  beaked  whales. 

OBJECTIVES 

To  achieve  the  long-term  goals,  a  number  of  scientific  and  technological  objectives  were 
identified.  These  included  the  following:  extending  existing  finite-element-method  (FEM)  code 
to  describe  acoustic  interactions  with  structures,  and  applying  this  to  a  virtual  beaked  whale  and 
mid-sonar  frequencies  in  the  range  1-10  kHz;  collecting  high-resolution  morphometric  data 
through  computerized  tomography  (CT)  scans  on  marine  mammal  specimens,  and  constructing 
finite-element  models  of  the  anatomy;  assigning  physical  properties  of  tissues;  benchmarking  the 
finite-element  code;  and  incorporating  the  extended  finite-element  code  and  morphometric  and 
physical-property  data  in  an  online  modeling  and  visualization  system  called  the  Virtual  Beaked 
Whale. 

APPROACH 

The  approach  followed  an  integrated  set  of  six  tasks,  which  are  briefly  described. 

Task  1 .  Development  of  a  finite-clement  method  to  model  acoustic  interactions:  This  task  was 
considered  as  a  structural-acoustic  problem,  in  which  most  tissue  groups  and  surrounding  water 
behave  as  acoustic  fluids  and,  bony  tissues  behave  as  elastic  solids.  A  beaked  whale  or  other 
marine  mammal  was  to  be  modeled  as  a  structure  represented  by  its  morphometry  (Task  2),  in 
which  each  anatomical  part  is  assigned  its  ow  n  set  of  physical  properties  (Task  3). 

Task  2.  Morphometry  and  meshing  the  three-dimensional  anatomy:  The  principal  source  of 
morphometric  data  were  CT  scans  performed  at  the  Woods  Hole  Oceanographic  Institution 
(WHOl)  Computerized  Scanning  and  Imaging  (CS1)  Facility.  Image  data  on  marine  mammal 
specimens  were  expressed  in  Digital  Imaging  and  Communications  in  Medicine  (DICOM) 
format.  Amira  (Visage  Imaging,  Inc.,  San  Diego,  CA)  visualization  software  was  used  for 
segmentation  and  surface  reconstruction.  The  automatic  mesh  generation  employed  tetrahedral 
elements. 

Task  3.  Physical  properties  of  tissues:  The  best  available  data  were  to  be  used  to  represent  the 
acoustically  important  properties  of  mass  density,  elastic  constants,  and  absorption  coefficients 
for  each  identified  internal  organ  or  other  body  part. 

Task  4.  Measuring  interactions  of  acoustic  fields  with  cetacean  specimens:  In  order  to  test  the 
FEM  code  (Task  1),  measurements  were  performed  of  the  internal  pressure  field  in  selected 
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tissues  in  an  instrumented  common  dolphin  cadaver,  also  referred  to  in  this  report  as  a  specimen. 
The  specimen  was  prepared  by  surgically  implanting  acoustic  sensors  in  the  form  of  tourmaline 
pressure  gauges;  CT-scanned  to  determine  sensor  location  and  morphometry;  then  acoustically 
measured  at  the  Naval  Surface  Warfare  Center  (NSWC)  Carderock  Division.  A  necropsy  was 
performed  within  hours  of  the  experiment. 

Task  5.  Testing  the  FEM  model:  Rigorous  testing  was  to  be  performed  by  comparison  with 
analytic  solutions  for  immersed  simple  objects.  These  solutions  were  developed  and  numerically 
realized  for  acoustically  absorptive  fluid  spheres  and  solid  elastic  spheres  in  a  lossy  immersion 
tluid.  Numerical  solutions  for  more  complicated  objects  were  also  identified. 

Task  6.  Virtual  Beaked  Whale:  This  interactive  online  modeling  and  visualization  system  was  to 
be  the  principal  deliverable  of  the  project.  It  was  to  incorporate  a  database  with  sets  of  whole- 
body  morphometric  data  (Task  2)  from  beaked  whales  and  other  species,  as  well  as  the 
respective  physical  properties  of  tissues  (Task  3).  It  would  also  allow  users  to  enter  other 
morphometric  and  physical-property  data  directly.  The  user  would  be  able  to  specify  an 
essentially  arbitrary  mid-frequency  sonar  signal.  The  output  would  consist  of  computed 
solutions  for  the  internal  pressure  and  displacement  fields  (Task  1)  at  user-specified  locations. 

WORK  COMPLETED 

Task  1.  Development  of  a  finite-element  method  to  model  acoustic  interactions 

2008:  Existing  finite-element-method  (FEM)  code  to  serve  as  the  core  of  the  modeling  system 
was  extended  in  several  ways.  (1)  A  perfectly  matched  layer  (Berenger  1094,  Turkel  and  Yefet 
1998)  was  implemented,  eliminating  extraneous  numerical  artifacts  of  wave  reflections  from  the 
boundary  of  the  computational  volume.  (2)  Frequency-domain  elastic  elements  were  realized  in 
code,  allowing  propagation  of  shear  waves  in  addition  to  ordinary,  longitudinal  acoustic  waves  in 
elastic  tissues  such  as  bone.  (3)  Interface  elements  were  added  to  transfer  loads,  namely  pressure 
and  stress,  between  acoustic  and  elastic  elements,  which  was  essential  for  representing 
interaction  effects  within  the  animal  model.  (4)  A  ten-node  tetrahedral  element  was 
implemented  in  the  code  allowing  for  quadratic  interpolation  of  the  basic  fields.  It  was  noted 
that  quadratic  elements  have  better  numerical  properties  than  linear  elements,  including  that  of 
numerical  dispersion,  or  wave-speed  effects  due  to  the  finite  size  of  elements,  as  well  as 
permitting  use  of  fewer  nodal  points  to  achieve  a  predefined  accuracy. 

2009:  An  essential  and  intricate  computational  problem  was  addressed.  This  concerned  the 
construction  of  software  capable  of  efficiently  solving  structural  acoustic  problems  with  millions 
of  degrees  of  freedom.  This  was  required  to  solve  problems  involving  acoustic  interactions  with 
heterogeneous  structures  and  generic  interfaces,  e.g.,  marine  mammals  in  their  natural  state  of 
immersion.  The  associated  FEM  code,  generieally  called  the  finite-element  solver,  was  now 
capable  of  solv  ing  general  structural-acoustic  problems.  This  incorporated  earlier  work 
implementing  a  perfectly  matching  layer,  realizing  acoustic  and  elastic  elements  and  interface 
elements  to  transfer  loads  between  these,  and  implementing  a  ten-node  tetrahedral  element  for 
quadratic  interpolation  of  the  basic  fields. 


5 


Task  2.  Morphometry  and  meshing  the  three-dimensional  anatomy 

2008:  Work  was  completed  on  several  operations  necessary  for  representing  the  animal  model  in 
a  three-dimensional  computational  space.  Anatomical  realizations  by  Amira  were  interfaced 
with  the  basic  meshing  software  and  FEM  software.  This  enabled  a  larger  computational  volume 
to  be  represented  where  the  perfectly  matching  layer  was  defined. 

The  meshing  process  is  illustrated  in  Fig.  1  using  morphometric  data  derived  from  CT  scans  on  a 
142-cm-long  common  dolphin  specimen.  Finite-element  models  based  on  tetrahedral  elements 
are  shown  in  Figs.  2-4.  The  scan  was  performed  at  8-mm  intervals,  hence  with  8-mm  slice 
thickness. 

A  seeond  specimen,  a  138-em-long  eommon  dolphin,  was  CT-scanned  at  the  WHOl  CS1  Facility 
at  much  higher  resolution  using  1-mm-thick  slices.  CT  data  on  a  living  bottlenose  dolphin  were 
received  from  the  U.S.  Navy  Marine  Mammal  Program. 

Images  of  tissues  were  also  identified  by  tissue  group  in  an  operation  called  labeling.  In  the  case 
of  the  same  142-cm-long  common  dolphin,  upper  respiratory  and  nearby  tissues  are  labeled  in 
Fig.  5.  Missing  lung  tissue,  due  to  the  resolution  of  the  8-mm  scan,  are  also  indicated. 

2009:  Another  essential  and  intricate  computational  problem  was  addressed.  This  concerned  the 
creation  of  finite-element  models  that  conform  to  tissue  interfaces.  This  required  the  use  of 
unstructured  meshes  as  opposed  to  rectangular,  Cartesian  grids.  CT  images  had  to  be  segmented; 
surfaces  between  different,  contiguous  tissues  had  to  be  triangulated;  and  the  encapsulated 
volumes  had  to  be  filled  with  volumetric  elements,  whieh  were  tetrahedral  for  the  sake  of 
computational  accuracy  and  efficiency.  These  volumetric  elements  were  geometrically 
constrained  to  avoid  extreme  acute  and  obtuse  shapes.  It  is  noted  that  commercially  available 
software,  sueh  as  Amira,  was  powerful  for  segmentation  and  surface  reconstruction,  but 
inadequate  for  automatic  generation  of  three-dimensional  meshes.  In  addition,  it  could  not 
generate  the  meshes  needed  for  defensible  acoustic  calculations,  such  as  the  creation  of  perfectly 
matched  layers  (PML),  among  other  things.  Using  customized  software,  three-dimensional 
visualizations  and/or  meshes  were  prepared  for  each  of  three  common  dolphins,  designated 
D_del40,  D_del42,  and  D  del52.  The  specimen  of  D  del52  was  used  in  the  experiment 
described  under  Task  4.  In  addition,  visualization  and  meshing  of  a  minke  whale  head  was 
commenced  in  anticipation  of  modeling  a  virtual  beaked  whale. 

2010:  The  principal  deliverable  of  the  project,  the  Virtual  Beaked  Whale  as  described  under  Task 
6,  required  a  morphometric  model  for  a  beaked  whale.  Potential  sources  of  such  data  included 
gross  anatomical  measurements  on  a  number  of  stranded  beaked  whales  and  computerized 
tomographic  (CT)  scans  of  the  heads  of  a  very  few,  different  beaked  whale  specimens.  Adult 
beaked  whales  are  simply  too  large  for  ordinary  CT-scanning.  To  remedy  this  shortcoming,  a 
synthetic  model  was  developed  based  on  the  following  sourees  of  data:  (i)  head  of  a  Cuv  ier's 
beaked  whale  ( Ziphius  cavirostris ),  a  sub-adult  male,  total  length  4  m,  mass  730  kg,  derived 
from  a  stranding  on  Cape  Cod  at  Sandy  Beach,  Cohasset,  MA,  on  4  April  2006,  with 
euthanization  by  the  Cape  Cod  Stranding  Network,  delivered  to  WHOl  for  dissection,  and  head 
CT-scanned  fresh  on  5  April  2006  at  WHOl,  and  (ii)  body  of  a  harbor  porpoise  ( Phocoena 


6 


phocoena),  an  adult  male,  length  105  cm,  mass  22. 1  kg,  dead  of  natural  causes  on  22  February’ 
2010,  found  in  Wellfleet,  MA,  delivered  to  WHOI  by  the  International  Fund  for  Animal  Welfare 
(1FAW),  and  scanned  fresh  the  following  day,  23  February  2010,  at  WHOI.  There  is  a 
presumption  that  the  anatomy  of  the  harbor  porpoise  body  resembles  that  of  a  Cuvier's  beaked 
whale,  at  least  to  within  a  scaling  factor.  7  he  corresponding  morphometric  data  were 
synthesized  by  reducing  the  head  of  the  beaked  whale  specimen  to  the  body  of  the  harbor 
porpoise,  rather  than  by  enlarging  the  body  of  the  harbor  porpoise  to  the  head  of  the  beaked 
whale.  More  particularly,  the  head  of  the  beaked  whale  was  scaled  and  stitched  to  the  body  of 
the  harbor  porpoise  within  the  Amira  environment. 

The  process  of  synthesizing  the  beaked  whale  model  is  illustrated  in  Figs.  6-8.  The  basic 
material  consists  of  CT  images  of  anatomical  cross  sections,  as  in  Fig.  6.  These  were  segmented, 
a  process  in  which  tissues  are  labeled,  allowing  automatic  linking  in  volumetric  rendering  across 
slices.  Ultimately  this  enables  sagittal  and  coronal  views  to  be  formed,  as  in  Fig.  7,  as  well  as 
transparent  views,  as  in  Fig.  8. 

Prior  to  modeling  the  sonar-induced  fields,  the  synthetic  beaked  whale  would  be  scaled  to  full 
beaked-whale  size.  In  the  present  ease  of  Cuvier's  beaked  whale,  the  overall  length  would  be  of 
order  4-8  m. 

201 1:  The  synthesis  process  was  completed  and  publication  was  initiated  on  the  CS1  website 
rhttn://csi.whoi.edu<'1. 

Task  3.  Physical  properties  of  tissues 

2000:  The  subject  of  physical  properties  of  marine  mammal  tissue  was  discussed  extensively.  It 
was  clear  that  current  knowledge  of  phy  sical  properties  was  mostly  derived  from  ex  situ 
measurements  rather  than  the  more  desirable  in  vivo  measurements.  It  was  also  appreciated  that 
the  physical  properties  of  some  tissues  are  influenced  by  pressure,  e.g.,  the  lungs  and  other  air- 
containing  tissues,  by  both  pressure  and  temperature,  as  in  the  case  of  the  melon. 

2010:  Physical  properties  of  beaked  whale  tissue:  This  was  a  continued  subject  of  discussion. 

The  basic  problem  was  the  lack  of  data,  further  exacerbated  by  the  use  of  frozen  and  thawed 
tissue  samples  for  the  few  measurements  that  have  been  undertaken  for  beaked  whale  tissue.  An 
unexpected  source  of  direct,  ex  situ  measurements  on  fresh  tissue  was  identified.  This  was  a 
series  of  measurements  of  the  longitudinal-wave  sound  speed  performed  by  D.  Chu  and  A. 
Lavery  on  tissue  samples  extracted  from  the  freshly  stranded  Cuvier's  beaked  whale  described 
above  under  Task  2.  Measurements  were  performed  with  a  unique  apparatus  (Chu  and  Wiebe 
2005),  which  had  earlier  been  applied  to  zooplankton  (Chu  et  al.  2000,  2003;  Wiebe  et  al.  2010). 
The  measured  beaked  whale  tissues  included  samples  from  the  ventral,  medial,  and  dorsal  parts 
of  the  melon;  blubber;  fat;  and  muscle.  The  Chu  and  Lavery  data  were  not  completely  analyzed 
follow  ing  their  collection  on  6  April  2006,  and  were  not  published. 

Gas  volume  conformational  changes  in  beaked  whales;  The  volume  of  gas  contained  within  air 
spaces  of  a  diving  and  ascending  odontocete  will  vary  depending  upon  depth,  functional  use,  and 
compliance  of  the  surrounding  tissues.  The  dynamic  relationship  between  gas  volumes  and 
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ambient  pressure  is  well  understood,  and  can  be  used  to  make  general  assumptions  about  the 
total  gas  volume  held  within  the  whale  at  any  depth  relative  to  the  total  volume  held  at  surface 
level.  However,  biological  material  is  not  uniform  and  therefore  variations  in  compliance  must 
be  accounted  for  when  examining  pressure-volume  relationships.  The  gas-containing  organs  can 
have  a  range  of  compliances  depending  upon  the  structure  of  the  interstitial  and  surrounding 
tissues.  Lungs,  for  example,  are  highly  elastic,  but  have  some  stretch  constraints  related  to  their 
internal  construction.  The  jointed  ribs  of  odontocetes  allow  additional  flexibility  for  lung 
collapse  compared  with  terrestrial  animals.  The  cartilage  surrounding  the  larynx  and  trachea, 
particularly  the  complete  circular  rings  found  in  whales,  makes  it  stiffer  than  the  lungs,  but  more 
pliant  than  the  skull.  Cranial  bony  sinuses  are  rigid  and  could  crack  under  pressure  changes  if 
the  ostia  allowing  gas  exchange  are  insufficient  or  blocked.  Odontocetes  have  thus  evolved  air 
sacs  outside  the  skull  that  are  free  of  the  constraints  of  rigid  bony  walls.  These  air  sacs  are 
probably  even  more  elastic  than  lungs,  particularly  as  they  have  highly  folded  external  contours 
and  no  tissues  crossing  through  their  lumen.  At  extreme  depths,  it  is  likely  that  these  spaces 
collapse  as  they  are  completely  evacuated  of  gas.  Some  species  of  deep  diving  beaked  whales 
have  evolved  reduced  or  absent  air  sacs,  probably  through  selection  against  these  non-functional 
spaces  at  full  collapse,  which  become  liability  for  infections.  There  are  several  major  functional 
conditions  under  which  odontocete  respiratory  spaces  may  undergo  conformational  changes 
during  a  dive:  breathing,  sound  production,  hearing,  and  valsalva.  It  is  clear  that  current 
know  ledge  of  beaked  whale  anatomy  is  inadequate  to  predict  gas  volume  conformational 
changes  with  depth  and  function. 

201 1 :  An  attempt  was  made  to  survey  the  literature  on  acoustically  significant  properties  of 
cetacean  tissues,  with  summary  in  Table  1.  These  properties  are  allowed  to  subsume  the 
conventional  physical  properties  of  mass  density,  sound  speed,  compressibility,  and  their 
pressure  and  temperature  dependences;  anatomical  properties,  including  morphology  and 
topology;  and  biochemistry,  including  molecular  composition.  The  source  literature  is  diverse, 
reflecting  the  many  interesting  aspects  of  cetacean  tissue  properties.  These  span  such  issues  as 
buoyancy  with  respect  to  temperature  and  pressure,  recognizing  that  lipids  undergo  a  phase 
change  within  the  ordinary  diving  range  of  sperm  whales  and  other  odontocetes;  role  of  blubber; 
audiology;  and  sound  generation,  including  echolocation.  The  importance  of  lipid  composition 
w  ith  tissue  location  is  appreciated  apropos  of  sound  propagation  (Koopman  et  al.  2003,  2006; 
Duggan  et  al.  2000).  Duggan  et  al.  (2009)  noted  the  relationship  of  carbon  chain  length  to  sound 
speed,  which  is  inverse  (Gouw  and  Vlugter  1967.  Hustad  et  al.  1971). 

It  is  observed  from  the  table  that  the  measurements  and  reported  data  are  quite  limited.  Notable 
exceptions  are  works  dating  from  Soldevilla  et  al.  (2005),  which  attempt  to  represent  the  physical 
properties  of  cetacean  tissues  in  a  comprehensive  manner  that  would  be  self-sufficient  for 
modeling  purposes.  In  Soldevilla  et  al.  (2005)  and  some  other  studies,  the  mass  density  of 
tissues  was  assumed  to  be  accessible  through  the  X-ray  absorption  measurements  that  are 
fundamental  to  computed  tomography  (CT),  and  conversion  from  Hounsfield  units  to  units  of 
mass  density  were  effected  through  calibrations  that  can  be  traced  to  measurements  on  phantom 
targets  of  known  mass  density.  The  represented  acoustic  properties  also  included  absorption. 

Determination  of  the  physical  properties  of  cetacean  tissues  represents  an  outstanding  problem 
that  has  attracted  attention  and  funding.  Thus  far,  tools  or  methods  to  accomplish  such 
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measurements  in  vivo  are  wanting.  Inference  of  tissue  properties  by  a  combined  program  of 
measurement  of  sonar-induced  fields  inside  instrumented  specimens  and  modeling  may  be  the 
most  effective  approach.  Certainly  the  issue  of  data  quality  is  important  when  justifying 
conclusions  being  drawn  from  modeling  studies  about  the  physical  effects  of  sonar. 

Task  4.  Measuring  interactions  of  acoustic  fields  w  ith  cetacean  specimens 

2008:  A  series  of  tests  were  performed  on  sensors  to  be  used  in  the  acoustic  experiments  with 
instrumented  marine  mammal  specimens  at  NSWC.  These  tests  established  that  tourmaline 
sensors  typically  used  at  shock-pressure  levels  can  also  be  used  at  sub-shock  levels,  and  that  new 
Lcnz-effect  transducers  can  generate  powerful  acoustic  signals  in  the  upper  part  of  the  mid¬ 
frequency  sonar  band.  Experience  was  gained  in  configuring  the  facility  for  experiments  with 
marine  mammal  specimens  through  measurements  performed  at  the  same  facility  on  a  30-kg  pig 
carcass.  This  carcass  was  instrumented  through  surgical  implantation  of  tourmaline  sensors,  and 
used  in  acoustic  tests  in  preparation  for  the  anticipated  marine  mammal  specimen  experiments. 

2009:  During  early  February  2009,  the  first  experiment  with  a  marine  mammal  specimen  was 
performed  at  NSWC  Cardcrock  Division.  The  specimen  was  that  of  a  young  adult  male  common 
dolphin  ( Delphinus  delphis),  1 19  kg,  21 1-cm  total  length,  that  had  stranded  and  died  on  a  Cape 
Cod  beach  two  months  previously.  In  the  interim  the  specimen,  designated  D  del52,  was  stored 
at  WHOl  in  a  freezer  at  -20°C.  CT  scans  were  obtained  in  the  frozen  state  and  following  thawing 
to  assess  overall  tissue  condition,  indicated  in  Fig.  9(a).  Tourmaline  pressure  gauges.  Fig.  9(b), 
were  surgically  implanted  in  the  nares,  left  lung,  right  lung,  external  right  lung,  melon,  along 
right  ear,  epaxial  muscle  near  the  dorsal  fin.  and  rectal  abdomen,  with  locations  summarized  in 
Table  11.  CT  scans  were  repeated,  Fig.  10,  to  determine  the  precise  anatomical  locations  of  the 
eight  gauges,  visualized  in  detail  in  Figs.  1 1  and  12.  The  specimen  was  transported  in  a  chilled 
state  to  NSWC  by  car. 

On  3  February  acoustic  calibration  data  were  collected  on  the  external  tourmaline  gauges 
mounted  on  the  suspension  frame  ensonified  with  a  pair  of  International  Transducer  Corporation 
(Santa  Barbara,  CA)  ITC-1032  transducers  arranged  in  a  dipole-like  configuration  with  a  15-cm 
rod  separating  the  two. 

On  4  February  the  specimen  was  suspended  within  a  PVC-pipe  frame.  Fig.  13.  and  immersed  in 
the  fresh  water  Test  Explosion  Pond,  Figs.  14-16,  with  specimen  at  2-m  depth.  It  was  exposed 
to  multiple  ensonifications  by  a  single  ITC-1032  spherical  transducer,  with  measured  source 
levels  described  in  Table  III.  This  was  excited  at  each  of  three  frequencies:  5,  7,  and  10  kHz,  by 
a  20-cycle  sinusoidal  burst  with  smooth  rise  and  fall  over  two  cycles,  for  each  of  three  source 
locations:  in  front  of  the  rostrum,  and  successively  on  the  right  and  left  sides  of  the  specimen, 
with  but  different  fore-aft  and  up-down  locations.  The  detailed  locations  arc  described  in  Table 
IV.  Data  were  recorded  simultaneously  on  15  channels  from  the  six  external  tourmaline  gauges, 
eight  implanted  gauges,  and  source  transducer,  over  a  total  of  5  s  for  each  frequency  and  source 
location,  with  4-Hz  pulse  repetition  frequency.  Each  signal  was  amplified  by  a  PCB  model 
422E01  inverting-charge  amplifier  and  digitized  at  100  kHz.  Overviews  of  the  acoustic  data  at  7 
kHz  are  presented  in  Fig.  1 7  for  the  source  on  the  right  side  of  the  specimen,  and  in  1  ig.  1 8  for 
the  source  in  front  of  the  specimen.  Details  from  Fig.  18  are  shown  in  Fig.  19  together  with  the 
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respective  results  of  band-pass  and  matched  filtering  of  the  several  signals.  Corresponding  data 
at  5  and  10  kHz  are  shown  in  Fig.  20  together  with  the  respective  processing  results.  Noise  data 
were  also  collected  passively,  without  source  excitation,  at  the  end  of  each  series. 

After  the  experiment,  the  specimen  was  removed  from  the  suspension  frame  and  transported  by 
car  to  the  Walter  Reed  Army  Institute  for  Research  (WRAIR),  where  a  necropsy  was  performed, 
led  by  Dr.  D.  R.  Ketten.  Tissues  were  observed  for  trauma  and  samples  collected  for  histological 
analysis.  A  WHOI  CSI-Scanning  and  Necropsy  Report  was  prepared. 

The  acoustic  data  were  comprehensively  reviewed  within  the  MATLAB  environment.  Noise 
was  prominent  at  60  Hz  and  at  sub-hertz  frequencies,  which  is  evident  in  the  total  signals 
presented  in  Figs.  17  and  18.  The  sub-hertz  noise  was  necessarily  non-stationary  owing  to  the 
short,  5-s  data  collection  periods.  Both  noise  types  were  removed  by  digital  band-pass  filtering 
with  a  bandwidth  equal  to  10%  of  the  center  frequency  of  the  transmit  signal.  The  effects  of 
such  filtering  supplemented  by  matched  filtering  are  shown  in  Fig.  21  for  the  case  of  the  source 
in  front  of  the  specimen  and  transmit  signal  centered  at  7  kHz.  Details  of  the  processing  are  also 
shown  in  Figs.  22  and  23.  Qualitative  conclusions  were  supplemented  by  a  detailed  quantitative 
analysis,  w'ith  results  presented  in  Table  V. 

2010:  In  the  experiment  performed  on  the  instrumented  specimen  of  a  young  adult  male  common 
dolphin  ( Delphinus  delphis)  at  NSWC  on  4  February  2009,  sonar-induced  internal  fields  were 
measured  with  tourmaline  gauges.  These  are  more  commonly  used  to  measure  high-amplitude, 
e.g.,  blast  fields  (Rogers  1997).  However,  they  are  evidently  suitable  for  measuring  low- 
amplitude,  sub-shock  fields  too.  During  the  experiment,  the  signals  at  the  output  of  the  gauges 
were  amplified  by  PCB  422  in-line  preamplifiers,  providing  a  voltage  gain  of  100.  Earlier 
performance  measurements  of  tourmaline  gauges  were  refined  and  extended.  In  particular,  two 
tourmaline  gauges  together  with  the  mentioned  preamplifiers  were  calibrated  using  a  calibrated 
Briiel  &  Kjaer  (Naerum.  Denmark)  B&K  8103  laboratory  hydrophone. 

2011:  Analysis  of  the  experiment  was  completed,  and  a  manuscript  was  prepared  and  submitted 
for  publication  (Foote  et  al.,  submitted). 

The  tourmaline  sensors  first  used  in  sub-shock  measurements  in  2008,  and  reported  above,  were 
documented  in  preliminary  form  (Morales  et  al.  2011). 

Task  5.  Testing  the  FEM  model 

2008:  Analy  tic  solutions  for  the  acoustic  interaction  of  plane  acoustic  w  aves  with  an  immersed, 
absorbing,  fluid  sphere  were  developed  and  rendered  in  code  for  use  in  validating  the  FEM  code. 
These  solutions  were  realized  numerically  (Foote  2007a,b)  for  a  50-mm-diameter  sphere 
ensonified  at  10  and  100  kHz.  The  immersion  medium  was  represented  as  water  of  mass  density 
1000  kg/m\  sound  speed  1500  m/s,  and  variable  absorption  in  the  range  0-10  dB/wavelength. 
The  sphere  was  represented  as  a  fluid  with  variable  mass  density'  over  the  range  500-2000  kg/m3 
to  achieve  density  contrasts  with  water  of  0.5-2,  variable  sound  speed  over  750-3000  m/s  to 
achieve  sound  speed  contrasts  with  water  of  0.5-2,  and  variable  absorption  coefficient  over  0-10 
dB/wavelength.  The  pressure  and  displacement  fields  were  computed  along  the  sphere  axis,  as 
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defined  by  direction  of  propagation  of  the  incident  wave;  transverse  to  this  axis  from  center;  and 
along  the  surface  from  the  forward  to  reverse  directions. 

2009:  The  FEM  code  described  in  Task  1  was  to  be  benchmarked  through  a  series  of  test 
computations.  (1)  Surface  fields  on  the  spherical  shape  due  to  plane-wave  ensonification  were 
computed  for  each  of  two  boundary  conditions:  pressure-release  and  rigid,  for  wavenumber- 
radius  products  over  the  range  0-15.  Farfield  backscattering  cross  sections  were  computed  on 
the  basis  of  the  surface  fields  and  compared  with  known  solutions.  (2)  Plane-wave-induced 
pressure  fields  in  a  fluid  sphere  with  relative  sound  speed  0.8  and  mass  density  equal  to  that  of 
the  immersion  medium  were  computed  for  wavenumber-radius  products  over  the  range  0-15. 

The  reflectivity  coefficient  was  computed  on  the  basis  of  the  surface  field  and  compared  with  the 
solution  due  initially  to  Anderson  (1950).  (3)  Internal  displacement  and  stress  fields  in  a  solid 
elastic  sphere  were  computed.  The  farfield  scattering  form  function  was  computed  as  a  function 
of  frequency  for  comparison  with  numerical  results  due  initially  to  Faran  (1951).  Relative  errors 
were  estimated  for  different  mesh  sizes  relative  to  the  acoustic  wavelength.  (4)  Internal 
displacement  and  stress  fields  induced  by  spherical  waves  incident  on  elastic  spherical  and 
cylindrical  shells  were  computed. 

Analytic  solutions  for  the  acoustic  interaction  of  both  plane  and  spherical  acoustic  waves  in  an 
absorbing  fluid  sphere  immersed  in  a  lossy  immersion  fluid  were  developed  and  rendered  in  code 
for  use  in  validating  the  FEM  code.  These  solutions  were  realized  numerically  for  a  50-mm- 
diameter  sphere  ensonified  at  10  and  100  kHz.  The  immersion  medium  was  represented  as  water 
of  mass  density  1000  kg/m3,  sound  speed  1500  m/s,  and  variable  absorption  in  the  range  0-10 
dB/wavelength.  The  sphere  was  represented  as  a  tluid  with  variable  mass  density  over  the  range 
500-2000  kg/m’  to  achieve  density  contrasts  with  water  of  0.5-2,  variable  sound  speed  over 
750-3000  m/s  to  achieve  sound  speed  contrasts  with  water  of  0.5-2,  and  variable  absorption 
coefficient  of  0  or  10  dB/wavelength.  The  pressure  and  displacement  fields  were  computed 
along  the  sphere  axis,  as  defined  by  direction  of  propagation  of  the  incident  wave;  transverse  to 
this  axis  from  center;  and  along  the  surface  from  the  forward  to  reverse  directions.  In  a  special 
additional  case,  that  of  a  focusing  sphere  with  fluid  mass  density  1900  kg/m  and  sound  speed 
857  m/s,  the  fields  were  also  computed  with  high  precision  at  selected  points  for  detailed 
comparison  and  assessment  of  computational  accuracy,  with  results  in  Fig.  24.  Preliminary 
documentation  is  contained  in  three  abstracts  (Foote  2007a, b;  Foote  and  Francis  2008). 

A  separate  solution  for  the  pressure  and  particle-velocity  fields  in  a  thin  spherical  shell  due  to  an 
external  point  source  is  presented  in  Fig.  25. 

2010:  Work  designed  for  benchmarking  computer  code  that  solves  the  full-diffraction  wave 
equation  inside  generally  irregular  and  heterogeneous  bodies  was  continued.  The  analytic 
solution  was  generalized  to  the  case  of  point  sources  in  the  nearfield  of  absorbing  fluid  spheres  in 
a  lossy  fluid  immersion  medium.  Detailed  results  suitable  for  benchmarking  are  presented  in 
Tables  VI  and  VII.  Tests  were  also  performed  using  boundary-element  and  finite-element  codes 
(Foote  et  al.  2009,  2010). 

Analytic  solutions  for  the  acoustic  interaction  of  both  plane  and  spherical  acoustic  waves  in  an 
absorbing  solid  elastic  sphere  immersed  in  a  lossy  immersion  tluid  were  also  developed  and 
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rendered  in  code.  These  solutions  were  realized  numerically  for  a  50-mm-diameter  sphere  of 
high-density  polyethylene  ensonified  by  a  plane  wave  of  amplitude  1  Pa  at  100  kHz.  The 
immersion  medium  was  represented  as  water  of  mass  density  1000  kg/m3,  sound  speed  1500  m/s, 
with  no  absorption  and  with  longitudinal-wave  absorption  0.4  dB/wavelength  and  transverse- 
wave  absorption  1 .2  dB/wavelength.  The  radial  stress  and  radial  displacement  fields  were 
computed  along  the  sphere  axis,  as  defined  by  direction  of  propagation  of  the  incident  wave, 
with  results  in  Fig.  26.  Preliminary  documentation  is  contained  in  an  abstract  (Foote  et  al.  2010). 

Task  6.  Virtual  Beaked  Whale 

2008:  Visualization  routines  were  developed  to  return  data  from  the  FEM  system  and  display 
these  on  a  webpage.  This  system  was  operational  in  a  pilot  mode  and  was  being  used  internally. 
Significant  achievements  included  specification  of  the  sonar  frequency  and  direction  of  the 
incident  plane  wave,  and  development  of  a  client-server  system  that  displayed  the  computational 
results  in  the  form  of  pressure  fields  in  user-defined  cross  sections  of  the  computational  domain. 
The  client-server  system  is  responsible  for  passing  data  among  the  different  components  of  the 
simulation  environment,  namely  meshing,  preprocessing,  and  FFM-analysis  programs.  The 
visualization  system  was  based  in  MATLAB,  and  was  capable  of  operating  with  meshes  of  up  to 
two  million  elements.  A  new,  successor  system  was  to  be  capable  of  operating  with  meshes  with 
more  than  20  million  elements. 

Work  was  commenced  on  developing  a  graphical  user  interface  for  the  Virtual  Beaked  Whale. 
Preliminary  results  are  shown  in  Fig.  27. 

2009:  A  user  group  composed  of  colleagues  from  the  international  bioacoustics  community  was 
defined  for  testing  the  Virtual  Beak  Whale  when  in  a  more  advanced  state. 

2010:  Modeled  results  for  the  pressure  magnitude  on  the  skeleton  of  the  common  dolphin 
specimen  D_del42  are  shown  in  Fig.  28. 

A  substantial  computing  platform  was  acquired  and  configured  to  support  the  Virtual  Beaked 
Whale.  This  was  an  HP  ProLiant  DL585  G6  server,  with  four  six-core  processors  and  four  500- 
GB  SATA  disks.  A  user’s  guide  for  operating  the  Virtual  Beaked  Whale  was  outlined. 

CONCLUSIONS 

Developing  an  interactive  online  modeling  and  visualization  system  to  predict  sonar-induced 
aeoustic  fields  inside  beaked  whales  as  well  as  inside  other  marine  mammals  is  a  formidable 
task.  It  was  the  overarching  goal  of  the  project  to  effect  this  system,  but  it  w  as  not  completed. 
However,  the  essential  composition  of  this  system  was  specified  in  detail,  recognizing  that  it 
would  be  necessary'  to  solve  several  tens  of  millions  of  equations  simultaneously. 

A  synthetic  model  of  Cuvier's  beaked  whale  was  developed  based  on  CT  scans  of  the  head  of  a 
freshly  dead  Cuvier’s  beaked  whale  and  CT  scans  of  the  body  of  a  harbor  porpoise.  These  data, 
in  segmented  form,  are  available  from  the  Woods  Hole  Oceanographic  Institution  (WHOl) 
Computerized  Scanning  and  Imaging  (CSI)  Facility  website  fhttp://esi. whoi.edu  1. 
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Source  literature  reporting  physical  property  values  of  cetacean  tissues  or  closely  related  topics 
has  been  summarized. 

Tourmaline  sensors  can  be  used  to  measure  relatively  weak,  sub-shock  acoustic  pressures. 

Sonar-induced  acoustic  pressure  fields  inside  a  common  dolphin  cadaver  have  been  measured. 
Both  amplitudes  and  times  of  flight  have  been  measured. 

Analytic  models  have  been  developed  for  benchmarking  FEM  code  based  on  the  interaction  of 
plane  harmonic  pressure  waves  with  lossy  fluid  and  solid  elastic  spheres. 

In  the  future,  experimental  measurements  of  sonar-induced  acoustic  fields  inside  instrumented 
marine  mammal  cadavers  may  be  useful  for  investigating  sound  transmission  pathways  and  for 
inferring  acoustic  properties  of  ensonfied  tissues. 
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Table  I.  Cetacean  tissue  properties  of  acoustic  significance. 


Tissue  Type 

Property  Investigated 

Approach 

Species  /  Name 

Nomenclature 

Khnow  ska  (1991) 

Nasal  tissues 

Phonic  fat  bodies 

Auditory  fat  bodies 

Blubber;  muscle 

Bony  tissues 

Anatomy 

Sound  speed 

Mass  density 

Pressure,  temperature 

dependences 

Biochemisty 

Methods  /  To  understand 

CT  *  computed  tomography 

GLC  *=  gas-liquid  chromatography 
MR1  =  magnetic  resonance  imaging 
TI  C  *  thin  layer  chromatography 

Source 

Phacoena  phocoena  / 
harbour  porpoise; 
Balaenaptera  musculus  / 
blue  whale;  Balaenoptera 
physalus  /  fin  whale 

X 

X 

Morphology  /  Thermal 
properties;  role  of  blubber 

Parry  (1949) 

Physeter  macracephalus ! 
sperm  whale 

X 

x 

X 

Density  analysis  /  Role  of 
spermaceti  organ  in 
buoyancy 

Clarke  (1970) 

Physeter  macracephalus  l 
sperm  whale 

X 

X 

TLC;  GLC;  lipid  &  fatty 
acid  analysis  / 1  ipid 
structure  of  spermaceti  organ 

Morris  (1973) 

Globicephala  melas  /  long- 
finned  pilot  whale 

X 

X 

X 

Dissection;  lipid  analysis; 
TLC;  GLC  /  Role  of  melon 
lipids  in  echolocation 

Wed  mid  et  al 
(1973) 

Tursiaps  t  run  cat  us  / 
bottlenose  dolphin 

X 

X 

Acoustic  time-of-  flight 
measurements  /  Sound 
transmission  in  head 

Norris  & 

Harvey 

(1974) 

Odontocetes 

X 

X 

Lipid  extraction;  TLC; 
in  fared  spectroscopy  /  Role 
of  melon  fats  in  echolocation 

Litchfield 

Greenberg 

(1974) 

Odontocetes 

X 

X 

X 

X 

X 

TLC;  lipid  &  fatty  acid 
analysis;  Wijs  method  / 
Echolocation  mechanisms 

I  itchilcld  ct 
til.  (1975) 

Physeter  macracephalus  l 
sperm  whale 

X 

X 

X 

Dissection;  T1  C:  G1  C  / 
Spermaceti  organ  function 

Morris  (1975) 

Porpoises  &  small  whales 

X 

X 

X 

X 

X 

X 

X 

Density;  T1  C;  lipid  analysis; 
ultrasound  /  Acoustic  lens 
for  echolocation 

Varanasi  ct 
al.  (1975) 

Balaenaptera  bare  alls  /  sci 
whale 

X 

X 

X 

Dissection;  histology;  lipid 
analysis  /  Sei  whale  blubber 

Ackman  et  al 
(1975) 

Glabicephala  melas  /  long- 
finncd  pilot  whale 

X 

X 

X 

Ultrasound  /  Role  of  melon 
in  echolocation. 

Blombcrg  & 
Jensen  (1976) 

Glabicephala  melas  /  long- 
finned  pilot  whale 

X 

X 

X 

TLC;  radioisotopes, 
ultrasound  /  Sound  velocities 
in  melon. 

Blombcrg  <fe 

Lindholm 

(1976) 

Kogia  b rev i ceps !  pygmy 
sperm  whale 

X 

X 

X 

Compositional  topography ; 
biochemistry;  GLC;  TLC  / 
Echolocation 

karol  et  al. 
(1978) 
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Table  I,  continued.  Cetacean  tissue  properties  of  acoustic  significance. 


Tissue  Type 

Property'  Investigated 

Approach 

Species  /  Name 

Nomenclature. 
Khnowska  (1991) 

Nasal  tissues 

Phonic  fat  bodies 

XA 

T3 

i 

<2 

c 

o 

t3 

< 

Blubber;  muscle 

Bony  tissues 

Anatomy 

Sound  speed 

Mass  density 

Pressure,  temperature 

dependences 

Biochemisty 

Methods  /  To  understand 

CT  "  computed  tomography 

GLC  =  gas-liquid  chromatographs 
MR1  =  magnetic  resonance  imaging 
TLC  =  thin  layer  chromatography 

Source 

Physeter  mocrocepholus  / 
sperm  whale 

X 

X 

X 

X 

X 

Dissection;  morphometry  / 
Structure  &  function  of 
spermaceti  organ 

Clarke 

(1978a) 

Physeter  mocrocephalus  / 
sperm  whale 

X 

X 

X 

X 

Head  oil  pressure  & 
temperature  measurements  / 
Buoy  ancy  function. 

Clarke 

(1978b) 

Physeter  macrocephalus  / 
sperm  whale 

X 

X 

X 

Modeling  /  Buoyancy 
control. 

Clarke 

(1978c) 

Physeter  mocrocepholus  / 
sperm  whale 

X 

X 

X 

X 

Dissection;  biochemisir}  of 
lipids  &  viscera  /  Functions 
of  fat  distribution 

l.ockyer 

(1991) 

Kogio  simus  /  dwarf  sperm 
whale. 

X 

X 

X 

Dissection;  veloeimctry  / 
Melon  function 

Goold  & 

Clarke  (2000) 

Kogio  breviceps  /  pygmy 
sperm  whale;  Kogio  simus  / 
dwarf  sperm  whale 

X 

X 

X 

X 

X 

Dissection,  thermal  anahsis 
/  Sound  production  & 
conduction;  buoyancy 
control 

Clarke  (2003) 

Odontocetes 

X 

X 

X 

X 

Dissection;  fatty  acid 
analysis:  GLC  /  Roles  of 
isovaleric  acid 

Koopinan  ct 
al.  (2003, 

2006) 

Ziphius  covirostris  f  Cuvier’s 
beaked  whale 

X 

X 

X 

X 

X 

X 

X 

X 

X 

Dissection;  CT  scans; 
veloeimctry  /  Physical 
property  values  determined 
for  modeling 

Soldcvilla  ct 
al.  (2005) 

Ziphius  covirostris  1  Cuv  ier’s 
beaked  whale 

X 

X 

X 

X 

X 

X 

X 

X 

CT  scans;  3D 

reconstruction;  veloeimctry ; 
finite  clement  modeling  / 

Phy  sical  effect  of  sound 

Krysl  el  al. 
(2006) 

Tursiops  truncotus  / 
bottlenose  dolphin 

X 

X 

X 

X 

X 

X 

X 

X 

CT  scans;  X-ray; 
veloeimctry  /  Postmortem 
changes;  tissue  structure; 
viscoelastic  properties 

McKenna  et 
al.  (2007) 

Ziphius  covirostris  /  Cuvier’s 
beaked  whale 

X 

X 

X 

X 

X 

X 

Dissection;  CT  scans  / 

Sound  transmission  & 
reception 

Cranford  ct  al. 
(2008b) 

Phocoena  phocoeno  / 
harbour  porpoise 

X 

X 

X 

Dissection;  CT  scans:  MR1 
Amira  /  Sound  generation; 
echolocation 

1  luggenberger 
ct  al.  (2009) 
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Table  II.  Tourmaline  gauge  placement  in  the  common  dolphin  specimen.  The  Cable  ID  denotes 
the  designation  of  the  gauge  signal  in  the  PCB  multi-channel  signal  conditioner  model  48 1 . 


Gauge 

No. 

Cable  ID 

Location 

Depth 

(mm) 

186 

M 

Nares 

120 

112 

H 

Left  lung 

130 

187 

N 

Right  lung 

110 

115 

1 

Thorax 

120 

90 

K 

Melon 

40 

188 

0 

Right  ear 

- 

171 

J 

Abdomen 

- 

05 

F 

Epaxial  muscle 

25 

Table  111.  Source  levels  of  the  1TC-1032  spherical  transducer  in  units  of  decibels  re  1  pPa 
at  1  m. 


Source 

Location 

5 

Frequency  (kHy.) 

7 

10 

1  Front 

158.6  ±0.3 

164.0  ±0.2 

169.6  ±0.1 

2  Right  side 

157.5  ±0.3 

162.1  ±0.1 

172.1  ±0.1 

3  Left  side 

155.3  ±0.2 

164.9  ±0.1 

168.9  ±0.1 

Table  IV.  Acoustic  source  locations  external  to  the  common  dolphin  specimen. 


Location 


Description 


1  10  cm  in  front  of  rostrum 

2  75  cm  aft  of  tip  of  rostrum,  67  cm  to  right  of 
specimen  centerline  as  viewed  from  above 

3  120  cm  aft  of  tip  of  rostrum,  95  cm  to  left  of 
specimen  centerline  as  viewed  from  above. 
35  cm  above  centerline 
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Table  V.  Statistical  characteristics  of  the  arrival  time.  Part  A,  and  amplitude.  Part  B,  of  the 
signals  in  five  tourmaline  pressure  gauges  used  as  receiving  sensors.  These  are  based  on  the 
arrival  time  and  amplitude  of  the  matched-filter  peak  for  the  direct-path  transmission  in  the 
external  receiving  sensor,  designated  Cable  B,  and  sonar-induced  pressure  fields  in  the  receiving 
sensors  implanted  in  the  epaxial  muscle.  Cable  F.  melon,  Cable  K,  nares.  Cable  M,  and  right 
lung,  Cable  N.  When  the  acoustic  source  was  on  the  left  side  of  the  specimen  and  transmitted  at 
5  kHz,  only  four  transmit  signals  were  recorded  on  Cable  Z,  hence  only  four  received  signals  in 
the  implanted  sensors  were  analyzed.  In  all  other  cases,  the  entire  set  of  20  transmit  and  received 
signals  were  analyzed. 


A.  Arrival  time  in  milliseconds 


Source  location 

Receiver  location 

5  kHz 

7  kHz 

10  kHz 

Front 

B:  External 

0.34  ±  0,01 

0.37  ±0.01 

0.37  ±0.00 

F:  Epaxial  muscle 

1.16  ±0.24 

2.79  ±  0.09 

2.06  ±  1.76 

K:  Melon 

0.23  ±  0.02 

0.25  ±0.01 

0.33  ±  0.00 

M:  Nares 

- 

0.42  ±  0.09 

0.62  ±  0.04 

N.  Right  lung 

- 

- 

Right  side 

B:  External 

0.26  ±0.01 

0.25  ±  0.00 

0  31  fO.OO 

F:  Epaxial  muscle 

0.58  ±0.02 

0.49  ±0.01 

0.57  ±  0.00 

K:  Melon 

0.80  ±  0.05 

0.59  ±  0.03 

0.67  ±  0.00 

M:  Nares 

0.60  ±  0.25 

0.73  ±0.05 

0.37  ±  0.03 

N:  Right  lung 

- 

0.84  ±0.1 5 

0.75  ±  0.04 

Left  side 

B:  External 

0.28  ±  0.01 

0.37  ±  0.00 

0.35  i  0.00 

F:  Epaxial  muscle 

0.86  ±0.10 

0.83  ±  0.04 

0.80  ±  0.02 

K:  Melon 

0.78  ±0.09 

1.08  ±0.03 

0.94  ±0.01 

M:  Nares 

- 

- 

1.04  t  0.19 

N:  Right  lung 

- 

- 

B.  Signal  magnitude  in  decibels 

Source  location 

Receiver  location 

5  kHz 

7  kHz 

10  kHz 

Front 

B:  External 

6.3  ±  0.3 

5.6  ±0.1 

5.6  ±0.1 

F:  Epaxial  muscle 

-20.3  ±  1.0 

-20.9  ±  0.9 

-25  7  £  0.8 

K:  Melon 

3.3  ±0.3 

3.9  ±  0.2 

8.2  ±  0.1 

M:  Nares 

- 

-20.1  ±0.7 

-20.3  ±  0.5 

N:  Right  lung 

- 

- 

Right  side 

B:  External 

8.6  ±  0.3 

8.9  ±0.1 

7.2  ±0.1 

F:  Epaxial  muscle 

-1.1  ±0.4 

-5.8  ±  0.2 

-3.8  ±0.1 

K:  Melon 

-3.2  ±0.3 

-7.3  ±  0.3 

-4.9  ±  0.2 

M:  Nares 

-15.5  ±0.8 

-18.0  ±0.8 

-20.4  ±  0.5 

N:  Right  lung 

- 

-20.5  ±  1 .0 

-21.7  ±0.4 

Left  side 

B:  External 

8.0  ±0.3 

5.6  ±  0.1 

6.1  ±  0.1 

F:  Epaxial  muscle 

-13.6  ±0.6 

-14.3  ±  0.3 

-13.3  ±  0.2 

K:  Melon 

-12.0  ±  1.2 

-10.5  ±0.3 

-7.6  ±0.2 

M:  Nares 

- 

- 

-28.5  ±  1.6 

N:  Right  lung 

- 

- 

- 
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Table  VI.  Internal  field  values  of  pressure  and  particle  displacement  due  to  point-source 
emonification  of  a  fluid  sphere  immersed  in  water  for  each  of  two  values  of  internal  absorption: 
(a)  ctiXi=0  dB,  and  (b)  aiXi=10  dB,  where  Xi  is  the  internal  wavelength.  The  sphere  diameter  2a 
is  50  mm,  mass  density  1900  kg/m\  and  sound  speed  857  m/s.  The  water  mass  density  is  1000 
kg/m'  and  sound  speed  1500  m/s.  The  point  source  is  located  75  mm  from  the  sphere  center.  Its 
source  strength,  qo  in  Eq.  (la),  is  1  Pa  m  and  its  harmonic  frequency  is  100  kHz.  The  field  point 
is  described  by  its  distance  r  from  the  sphere  center  and  polar  angle  0,  which  is  the  angle  between 
the  vector  to  the  field  point  and  longitudinal  axis  of  the  sphere.  In  the  implied  coordinate 
system,  the  source  is  located  at  ro=75  mm  and  0=180  deg. 


(a)  a\!  X\  -  0  dB/m 


r/a 

0 (deg) 

Pin.  (Pa) 

ur  (pm) 

U(|  (pm) 

Real 

lmag 

Mag 

Real 

lmag 

Mag 

Real 

lmag 

Mag 

-1.0 

180 

-12.2831 

17.4343 

21.3267 

6.8180 

7.4401 

10.0916 

0 

0 

0 

-0.5 

180 

17.4377 

-24.7762 

30.2974 

-13.3784 

-6.8501 

15.0301 

0 

0 

0 

0.0 

- 

1.1625 

25.2769 

25.3036 

-22.8015 

0.3063 

22.8036 

0 

0 

0 

0.5 

0 

3.5059 

-27.7566 

27.9771 

22.7649 

-15.4230 

27.4975 

0 

0 

0 

1.0 

0 

92.4330 

84.0886 

124.9590 

-59.6620 

20.0940 

62.9549 

0 

0 

0 

0.5 

30 

-35.6583 

-1.7023 

35.6989 

-5.8754 

-23.3146 

24.0435 

1.4424 

13.1947 

13.2733 

0.5 

60 

16.4382 

-14.0547 

21.6275 

2.8506 

6.7643 

7.3405 

-21.6775 

-20.3666 

29.7441 

0.5 

90 

14.9343 

19.0374 

24.1962 

2.0098 

- 1 .9264 

2.7840 

13.2588 

-1  1.8924 

17.8108 

0.5 

120 

-21.0330 

5.6706 

21.7840 

3  8198 

11.9976 

12.5910 

2.5568 

16.8144 

1 7.0077 

0.5 

150 

21.5088 

-3.7667 

21.8361 

-3.2086 

-19.0306 

19.2991 

-5.4500 

-8  8940 

10.4310 

fb)  a,/X,  -  10  dB/m 

Pin.  (Pa) 

ur  (pm) 

ue  (pm) 

r/a 

0 (deg) 

Real 

lmag 

Mag 

Real 

lmag 

Mag 

Real 

lmag 

Mag 

-1.0 

180 

-9.1557 

18.4622 

20.6077 

16.5499 

12.1381 

20.5239 

0 

0 

0 

-0.5 

180 

1.0245 

-4.0501 

4.1777 

-3.8124 

-1.6951 

4.1723 

0 

0 

0 

0.0 

- 

-0.0625 

0.8516 

0.8539 

-0.8321 

-0.2065 

0.8573 

0 

0 

0 

0.5 

0 

-0.0861 

-0.3580 

0.3682 

0.1392 

-0.1555 

0.2087 

0 

0 

0 

1.0 

0 

-0.2023 

1.1073 

1.1256 

0.3110 

0.3398 

0.4606 

0 

0 

0 

0.5 

30 

-0.2370 

-0.1888 

0.3030 

0.1210 

-0.3134 

0.3360 

-0.2170 

0.0531 

0.2234 

0.5 

60 

0.8397 

-0.4586 

0.9568 

0.2491 

0.1052 

0.2704 

-0.3270 

-0.8924 

0.9504 

0.5 

90 

0.9890 

1 .3942 

1.7094 

0.4617 

-0.0396 

0.4634 

1 .4422 

-0.8068 

1.6525 

0.5 

120 

-2.6234 

0.6554 

2.7040 

0.0083 

1.6722 

1 .6722 

0.2995 

2.1047 

2.1259 

0.5 

150 

3.6525 

-0.7000 

3.7190 

-0.0197 

-3.3134 

3.3134 

-0.1461 

1.6761 

1 .6824 
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Table  VII.  Internal  field  values  of  pressure  and  particle  displacement  due  to  plane-wave 
emonificalion  of  a  fluid  sphere  immersed  in  water  for  each  of  two  values  of  internal  absorption: 
(a)  ai^i=0  dB,  and  (b)  ai^i=10  dB.  where  Xi  is  the  internal  wavelength.  The  sphere  diameter  2a 
is  50  mm,  mass  density  1900  kg/m\  and  sound  speed  857  m/s.  The  water  mass  density  is  1000 
kg/m'  and  sound  speed  1500  m/s.  The  plane  wave  propagates  along  the  longitudinal  axis  of  the 
sphere.  Its  amplitude,  po  in  Eq.  (5a),  is  1  Pa  and  its  harmonic  frequency  is  100  kHz.  The  field 
point  is  described  by  its  distance  r  from  the  sphere  center  and  polar  angle  0,  which  is  the  angle 
between  the  vector  to  the  field  point  and  longitudinal  axis  of  the  sphere.  In  the  implied 
coordinate  system,  the  source  is  located  at  infinity  with  0=180  deg. 


(a)  a,/  X)  =  0  dB/m 


r/a 

0 (deg) 

Pirn  (Pa) 

ur  (pm) 

Uu  (pm) 

Real 

lmag 

Mag 

Real 

lmag 

Mag 

Real 

lmag 

Mag 

-1.0 

180 

-0.9584 

2.6523 

2.8202 

1.5752 

-0.0441 

1.5759 

0 

0 

0 

-0.5 

180 

-0.1511 

-2.0656 

2.0711 

-1.5746 

-0.7560 

1.7467 

0 

0 

0 

0.0 

- 

0.0872 

1.8958 

1.8978 

-1.7077 

0.0773 

1.7094 

0 

0 

0 

0.5 

0 

0.2950 

-3.0623 

3.0765 

3.0427 

-0.6949 

3.1210 

0 

0 

0 

1.0 

0 

6.2740 

-5.1387 

8.1098 

1.7021 

6.2059 

6.4351 

0 

0 

0 

0.5 

30 

-2.2938 

2.1729 

3.1596 

-1.1312 

-0.9730 

1.4921 

1 .4993 

1 .3596 

2.0240 

0.5 

60 

-0.4380 

-1.4250 

1 .4908 

0.1945 

0.1891 

0.2713 

-1.7900 

1  0743 

2.0876 

0.5 

00 

1.6826 

-0.5235 

1.7622 

-0.0992 

-0.5083 

0.5179 

-0.1970 

-0.9603 

0.9803 

0.5 

120 

-0.7858 

1.1184 

1.3669 

0.8486 

0.4313 

0.9519 

0.9485 

0.5828 

1.1 132 

0.5 

150 

1.1298 

-0.5197 

1.2436 

-0.5815 

-1.1602 

1.2978 

-0.0617 

-0.2152 

0.2238 

(b)  a,/  X,  = 

10  dB/m 

PiM  (Pa) 

ur  (pm) 

Ufl  (pm) 

r/a 

0 (deg) 

Real 

lmag 

Mag 

Real 

lmag 

Mag 

Real 

lmag 

Mag 

-1.0 

180 

-0.4436 

0.9290 

1 .0295 

0.8426 

0.5797 

1 .0228 

0 

0 

0 

-0.5 

180 

0.0570 

-0.2399 

0.2465 

-0.2276 

-0.0917 

0.2454 

0 

0 

0 

0.0 

- 

-0.0047 

0.0639 

0.0640 

-0.0628 

-0.0135 

0.0643 

0 

0 

0 

0.5 

0 

-0.0156 

0.0172 

0.0232 

-0.0034 

0.0307 

0.0309 

0 

0 

0 

1.0 

0 

0.1822 

-0.0517 

0.1894 

0.0008 

-0.0997 

0.0997 

0 

0 

0 

0.5 

30 

-0.0587 

0.0159 

0.0608 

-0.0249 

-0.0235 

0.0342 

-0.0005 

0.0570 

0.0570 

0.5 

60 

-0.0309 

-0.0998 

0.1045 

-0.0101 

-0.0144 

0.0176 

-0.1025 

0.0224 

0.1049 

0.5 

90 

0.1459 

-0.0382 

0.1508 

0.0050 

-0.0682 

0.0684 

-0.0248 

-0.1322 

0.1345 

0.5 

120 

-0.1076 

0.1660 

0.1978 

0.1038 

0.1022 

0.1456 

0.1069 

0.0799 

0.1334 

0.5 

150 

0.2119 

-0.0977 

0.2334 

-0.0579 

-0.2085 

0.2164 

-0.0308 

-0.0793 

0.0851 
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Fig.  1.  Construction  of  the  finite-element  mesh  for  a  142-cm-long  common  dolphin .  From  CT 
images  (a)y  surface  reconstructions  of  tissues  can  he  obtained  (b).  These  surfaces  are  sub- 
sampled  and  triangulations  of  the  interfaces  are  created.  From  these  triangnlationsy  a 
tetrahedral  mesh  is  created  (c)y  illustrated  in  greater  detail  in  Figs.  2-4.  An  external  meshy 
which  is  necessary  for  acoustic  computations ,  is  constructed  and  joined  with  the  previous 
mesh.  An  outline  of  this  is  shown  in  (d)y  with  different  perfectly  matched  layer  blocks  shown 
in  color.  The  meshed  dolphin  within  the  computational  volume  is  shown  in  (e).  A  detail  of 
the  triangulation  in  the  surface  of  the  external  mesh  is  shown  in  (f). 

Credit:  G.  Feijooy  D.  Ketten. 
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Fig.  2 .  Finite-element  representation  with  tetraheda  of  a  142-cm-long  common  dolphin  and 
external ,  computational  volume ;  ///  two  views .  foj  The  lungs  are  indicated  in  blue,  (b)  The 
skeleton  is  indicated  in  gray  and  the  lungs ,  in  green .  Credit:  G.  Feijoo ,  Z).  Ketten. 
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(b) 


tig.  3.  Details  of  tetrahedral  meshing  for  the  common  dolphin  specimen  in  Fig.  2. 

Credit:  G.  Feijoo,  D.  Ketten. 
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(b) 

Fig-  4.  Further  details  of  tetrahedral  meshing  for  the  common  dolphin  specimen  in  Figs.  2 

and  3.  Credit:  G.  Feijoo,  D.  Ketten. 
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Light  blue  area 
indicates 
approximate 
anatomical 
outline  of  missing 
lung  region 


Higher  CT  scan  resolution  will  provide  better  discrimination 
between  thin  layers  of  air  filled  lung  and  adjoining  soft  tissues 


Fig.  5.  Three-dimensional  surface  rendering,  based  on  CT  scans  visualized  with  Amira,  of  the 
lower  respiratory  tract  in  a  142-cnt-long  common  dolphin  specimen,  with  superimposed 
transverse  CT  section.  Upper  panel:  lower  respiratory'  tract  and  nearby  tissues  labeled.  Lower 
panel:  missing  lung  volumes  indicated.  Credit:  J.  Reidenberg. 
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Fig.  6.  Transaxial  two-dimensional  slice  through  the  thoracic  region  of  the  synthetic  model  of 
a  beaked  whale,  with  lungs  outlined  in  red,  and  skeleton,  consisting  of  ribs,  one  vertebrate, 
and  bones  in  the  pectoral  flippers,  outlined  in  green.  Credit:  D.  Ketten,  S.  Cramer. 
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Fig.  7.  Axial,  sagittal,  and  corona!  views  of  the  synthetic  model  of  a  beaked  whale,  together 
with  a  dorsal  transparent  view  highlighting  the  skeleton  and  air  passages. 

Credit:  D.  Ketten,  S.  Cramer. 
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Fig.  8.  Left  lateral  transparent  view  of  the  synthetic  model  of  a  beaked  whale,  highlighting  the 
skeleton,  air  passages,  and  exterior  surface.  Credit:  D.  Ketten,  S.  Cramer. 
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Fig.  9.  (a)  The  experimental  subject,  a  1 19-kg,  21 1-ctn-long  common  dolphin  specimen, 
following  implantation  of  several  tourmaline  pressure  sensors  (b)  at  the  WHOI  CSl  Facility. 

Credit:  D.  Ketten,  S.  Cramer,  J.  Aruda. 
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(b) 


Fig.  10.  The  experimental  subject  during  implantation  of  tourmaline  pressure  sensors  at  the 
WHOI CS1  Facility,  (a)  with  use  of  the  CT  scanner  to  confirm  sensor  placements,  and 
(b)  following  completion  of  sensor  placements.  Credit:  D.  Ketten,  S.  Cramer ,  J.  Aruda. 
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Fig.  1 1.  Gauge  placements  in  the  common  dolphin  specimen  according  to  a  virtual  ray  tracing 
(VRT)  three-dimensional  image,  which  also  includes  the  lungs ,  skeleton,  and  skin  for 
reference.  Credit:  D.  Ketten,  S.  Cramer,  J.  Aruda. 
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Fig.  12.  Head  gauge  placements  in  the  melon,  nares,  and  ear  of  the  common  dolphin 
specimen  whose  virtual  ray  tracing  (VRT)  three-dimensional  image  is  shown  in  Fig.  II,  which 
also  includes  the  lungs,  skeleton,  aud  skin  for  reference. 

Credit:  D.  Ketten,  S.  Cramer,  J.  Arudu. 
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Fig.  13.  The  experimental  subject,  a  119-kg,  21 1-cm-long  common  dolphin  specimen, 
following  implantation  of  tourmaline  sensors  at  the  WH01 CSI  Facility,  suspended  from  a 
PVC-pipe  frame  at  NSWC  Carderock  Division  on  4  February  2009.  Credit :  D.  Ketten,  S. 
Cramer,  K.  Rye,  W.  11.  Lewis,  M.  Hastings. 


Fig.  14.  Explosion  Test  Pond  facility  at  NSWC  Carderock  Division:  bird’s-eye  view. 
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Fig.  IS.  Specimen  in  its  suspension  frame  being  lowered  into  the  pond  by  a  crane.  Credit:  D. 
Ketten,  S.  Cramer,  K.  Rye,  W.  H.  Lewis,  M.  Hastings. 
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Fig.  16.  The  experimental  subject,  a  119-kg ,  21  Fern-long  common  dolphin  specimen , 
following  implantation  of  tourmaline  sensors  at  the  WHOI CSI  Facility ,  suspended  from  a 
PVC-pipe  frame  at  NS WC  Carderock  Division  on  4  February  2009 ,  and  immersed  in  the 
NSWC  Explosion  Test  Pond  facility  for  acoustic  measurement ,  lower  frame. 

Credit:  D.  Ketten ,  S.  Cramer ;  A".  /?>><?,  W7.  //.  Lewis,  M.  Hastings. 
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Fig .  17.  Overview  of  acoustic  data  collected  during  the  experiment  with  the  common  dolphin 
specimen  at  NSWC  Carderock  Division  on  4  February  2009  with  the  ITC-1032  spherical 
transducer  on  the  right  side  of  the  specimen  and  transmit  signal  centered  at  7  kHz.  The  total 
data  collection  time  is  5  s;  the  pulse  repetition  frequency  is  4  Hz.  The  total  signal ,  including 
background  noise ,  is  shown  in  gray.  The  signal  resulting  from  digital  bandpass  filtering  with 
bandwidth  equal  to  one-third  of  the  center  frequency  is  shown  in  red.  Cables  A ,  B,  C,  D,  E, 
and  G  are  connected  to  external  sensors.  Implanted  sensors  are  connected  to  Cable  F  for  the 
epaxial  muscle  near  the  dorsal  fin ,  Cable  H  for  the  left  lung ,  Cable  1  for  the  thorax ,  Cable  J 
for  the  abdomen ,  Cable  Kfor  the  melon,  Cable  M for  the  nares,  Cable  N  for  the  right  lung, 
and  Cable  O  for  the  right  ear.  Credit:  Y.-T.  Lin,  M.  Hastings,  K.  Rye ,  K.  Foote. 
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Fig.  18.  Overview  of  acoustic  data,  as  in  Fig.  17,  but  with  the  spherical  transducer  in  front  of 
the  specimen,  again  with  the  transmit  signal  centered  at  7  kHz. 
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Fig.  19.  Acoustic  data  from  Fig.  18  together  with  results  of  hand-pass  and  matched-  filtering 
for  external  sensor  B  in  the  line  of  sight  of  the  spherical  transducer  in  front  of  the  specimen 
and  internal  sensor  K  implanted  in  the  melon,  with  the  transmit  signal  centered  at  7  kHz- 
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Fig.  20.  Details  of  acoustic  data  for  externa!  sensor  B  in  the  line  of  sight  of  the  spherical 
transducer  in  front  of  the  specimen  and  internal  sensor  K  implanted  in  the  melon,  with  the 
transmit  signal  centered  at  (a)  5  kHz  and  (b)  10  kHz. 
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Fig.  21.  Overview  of  acoustic  data  in  Fig .  18,  with  the  spherical  transducer  in  front  of  the 
specimen  and  transmit  signal  centered  at  7  kHz »  but  following  band-pass  and  matched-filter 

processing . 


45 


Transmitted 

signal 


I 


I 


Time  (ms) 


Fig.  22.  Detail  of  processed  signals  derived  from  Fig.  17,  with  the  spherical  transducer  on  the 
right  side  of  the  specimen  and  the  transmit  frequency  centered  at  7  kHz*  The  electronic 
transmit  signal,  shown  in  the  upper  frame,  was  recorded  on  the  channel  marked  Cable  Z. 
Corresponding  signals  received  in  the  epaxial  muscle,  Cable  F,  and  melon.  Cable  K,  are 
shown  with  their  7-kHz  center  frequency  in  the  middle  and  lower  frames ,  respectively.  The 
simple  functions  shown  in  the  middle  and  lower  frames  are  the  respective  matched-filter 
outputs,  derived  according  to  the  method  in  Task  4. 


Amplitude 
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Fig.  23.  Matched-filter  outputs  of  5-  and  7-kHz  signals  received  in  the  epaxial  muscle  and 
melon  when  the  acoustic  source  was  positioned  in  front  of  the  specimen.  Each  displayed 
matched-filter  output  is  the  average  of  the  20  corresponding  outputs. 
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Fig.  24.  Pressure  and  longitudinal  particle  displacement  amplitudes  modeled  along  the 
longitudinal  axis  of  an  immersed  50-mm-diameter  fluid  sphere  at  100  kHz  due  to 
ensonification  by  a  spherical  wave  emanating  from  a  point  source  of  strength  1  Pa  m  at  z=-75 
mm  on  the  same  axis  (a,c),  and  a  plane  wave  of  amplitude  1  Pa  propagating  in  the  positive  z- 
direction  (b,d).  The  mass  density  of  the  fluid  sphere  is  1900  kg/m3,  sound  speed,  857  m/s,  and 
absorption  coefficient,  «/,  0  or  1  dB  per  internal  wavelength.  The  mass  density  of  the 
immersion  medium  is  1000  kg/m3,  sound  speed,  1500  m/s,  and  absorption  coefficient,  0  dB/m. 
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Fig.  25.  Pressure  field  in  water  (colors)  and  velocity  field  (sub-sampled  and  displayed  using 
arrows)  modeled  in  a  spherical  shell  due  to  a  monopole  source.  The  thin  spherical  shell  was 
not  rendered  to  permit  visualization  of  the  velocity  field.  Credit:  G.  Feijoo. 
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Fig.  26.  Radial  stress  (top)  and  displacement  (bottom)  amplitudes  modeled  inside  a  50-mm- 
diameter  sphere  of  higli-density  polyethylene  due  to  an  incident  100-kHz  plane  wave  of 
amplitude  1  Pa.  The  mass  density  and  longitudinal-  and  transverse-wave  sound  speeds  of  the 
HDPE  are  957  kg/m3,  2450  m/s,  and  950  m/s,  respectively.  The  immersion  medium  is  water, 
with  mass  density  1000  kg/m3  and  sound  speed  1500  m/s.  The  computations  have  been 
repeated  assuming  no  losses  in  the  sphere  and  assuming  the  following  values  for  the  product 
of  absorption  coefficient  and  wavelength:  0.40  and  1.20  dB  for  the  respective  longitudinal- 
and  transverse-wave  sound  speed.  Credit:  D.  T.  I.  Francis. 
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Fig.  27.  Virtual  Beaked  Whale:  development  of  a  graphical  user  interface ,  with  field  values  of 
the  function  f(x,y9z)  ~xT-y2-t  shown  on  the  computational  boundaries  in  three  orthogonal 
planes  and  on  (a)  two  internal \  orthogonal  slices ,  and  (b)  an  iso-surface. 

Credit:  G.  Feijoo 9  L.  Cohen. 


(C) 


Fig.  28.  Virtual  Beaked  Whale:  initial  simulations  using  D_del42  specimen  scanned  in  2006. 
The  lungs  were  not  modeled,  but  assumed  to  have  the  same  properties  as  the  surrounding 
water.  The  incident  sonar  signal  is  a  5-kHz  plane  wave  moving  in  the  negative  z-direction. 
The  magnitude  of  pressure  is  simulated  on  (a)  the  skeleton,  (b)  the  skeleton  and  in  a  saggital 
cut,  and  (c)  the  sagittal  cut  only.  Credit:  G.  Feijoo. 


